Conjugate points in Euler's elastic problem* 



Yu. L. Sachkov 
Program Systems Institute 
Russian Academy of Sciences 
Pereslavl-Zalessky 152020 Russia 
E-mail: sachkov@sys.botik.ru 

February 5, 2008 



Abstract 

For the classical Euler's elastic problem, conjugate points are de- 
scribed. Inflectional elasticae admit the first conjugate point between 
the first and the third inflection points. All the rest elasticae do not have 
conjugate points. 

Keywords: Euler elastica, optimal control, conjugate points, exponential 

mapping 

Mathematics Subject Classification: 49J15, 93B29, 93C10, 74B20, 
74K10, 65D07 



* Work supported by the Russian Foundation for Basic Research, project No. 05-01-00703-a. 



1 



Contents 



1 Introduction 3 

2 Conjugate points, Morse index, and MeisIov index 4 

2.1 Optimal control problem and Hamiltonians 4 

2.2 Second variation and its Morse index 5 

2.3 Exponential mapping 7 

2.4 Maslov index of a curve in Lagrange Grassmanian 7 

2.5 Morse index and Maslov index 8 

2.6 Preliminary remarks on Euler's problem 10 

3 Conjugate points on inflectional elasticae 11 

3.1 Preliminary lemmas 12 

3.1.1 Roots of the function oq 12 

3.1.2 Roots of the function ao + ai + 02 15 

3.2 Bounds of conjugate time 17 

3.3 Upper bound on cut time 24 

4 Conjugate points on non- inflectional elcisticae 25 

4.1 Preliminary lemmas 25 

4.2 Absence of conjugate points on non-inflectional elasticae 27 

4.3 Absence of conjugate points for special cases 28 

5 Final remarks 28 
List of Figures 30 
References 31 



2 



1 Introduction 



This work is devoted to the study of the following problem considered by Leon- 
hard Euler [7]. Given an elastic rod in the plane with fixed endpoints and 
tangents at the endpoints, one should determine possible profiles of the rod 
under the given boundary conditions. Euler's problem can be stated as the 
following optimal control problem: 

a; = cos 6*, (1.1) 
y = sinO, (1.2) 
= u, (1.3) 

q^{x,y,e) e M = M.lyX S^, ueR, (1.4) 
q{0) = qo = {xQ,yo,Oo), q{ti) = qi = {xi,yi,9i), fixed, (1.5) 

1 /"'^ 

J=- u^{t)dt^ min, (1.6) 

2 Jo 

where the integral J evaluates the elastic energy of the rod. 

This paper is an immediate continuation of our previous work [11], which 
contained the following material: history of the problem, description of attain- 
able set, proof of existence and boundedness of optimal controls, parametriza- 
tion of extremals by Jacobi's functions, description of discrete symmetries and 
the corresponding Maxwell points. In this work we widely use the notation, 
definitions, and results of work [11]. 

Leonhard Euler described extremal trajectories of problem (1.1) (1.6), their 
projections to the plane {x, y) being called Euler elasticae. Although, the ques- 
tion of optimality of elasticae remained open. Our aim is to characterize global 
and local optimality of Euler elasticae. Short segments of elasticae are optimal. 
The main result of the previous work [11] in this direction was an upper bound 
on cut points, i.e., points where elasticae lose their global optimality. In this 
work we describe conjugate points along elasticae; we obtain precise bounds for 
the first conjugate point, where the elasticae lose their local optimality. 

Each inflectional elastica contains an infinite number of conjugate points. 
The first conjugate point occurs between Maxwell points; visually, the first 
conjugate point is located between the first and the third inflection point of the 
elastica. 

All the rest elasticae do not contain conjugate points. 

Notice that Max Born proved in his thesis [5] that if an elastic arc is free of 
inflection points, then it does not contain conjugate points, so in this part we 
repeated Max Born's result. Although, our method of proof is more flexible, 
and we believe that it will be useful for the study of conjugate points in other 
optimal control problems. 

This work has the following structure. In Sec. 2 we review some basic facts 
of the theory of conjugate points along regular extremals of optimal control 
problems. These facts are rather well known, but are scattered through the 
literature. The main facts of this theory necessary for us are the following: 
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(1) An instant t > is a conjugate point iff the exponential mapping for time t 
is degenerate; (2) Morse index of the second variation of the endpoint mapping 
along an extremal is equal to the number of conjugate points with account of 
their multiplicity: (3) Morse index is equal to Maslov index of the curve in a 
Lagrange Grassmanian obtained by linearization of the flow of the Hamiltonian 
system of Pontryagin Max;imum Principle; (4) Maslov index is invariant under 
homotopics of extremals provided that their endpoints arc not conjugate. We 
apply this theory for description of conjugate points in Euler's problem. In Sec- 
tion 3 we obtain estimates for the first conjugate point on inflectional elasticae. 
Moreover, we improve our result of work [11] on the upper bound of cut time 
on inflectional elasticae. In Section 4 we show that all the rest elasticae do not 
contain conjugate points. In Section 5 we present some final remarks on results 
obtained in this paper and discuss their possible consequences for future work. 

In this work we use extensively Jacobi's functions, see [8, 13]. We apply 
the system "Mathematica" [14] to carry out complicated calculations and to 
produce illustrations. 

2 Conjugate points, Morse index, 
and Maslov index 

In this section we recall some basic facts from the theory of conjugate points in 
optimal control problems, see [1-4,12]. 

2.1 Optimal control problem and Hamiltonians 

We consider an optimal control problem of the form 



where M is a finite-dimensional analytic manifold, f{q,u) and (p{q,u) are an- 
alytic in {q, u) families of vector fields and functions on M depending on the 

control parameter u E U, and U an open subset of M™. Admissible controls axe 
u{-) G Loo[0,ti], and admissible trajectories q{-) are Lipschitzian. Let 

hu{X) = {X,f{q,u))-ip{q,u), \eT*M, q = n{\)eM, u e U, 

be the normal Hamiltonian of PMP for the problem (2.1)-(2.3). Fix a triple 
{u{t), Xt,q{t)) consisting of a normal extremal control u{t), the corresponding 
extremal A/, and the extremal trajec;tory q{f.) for the problem (2.1)-(2.3). 
In the sequel we suppose that the following hypothesis holds: 

(HI) For all A e T*M and u & U , the quadratic form ^ ^ (A) is negative 
definite. 



Q = f{Q,u) 
9(0) = qo, 



qGM, u e C/ C M", 
q{ti) = qi, ti fixed, 



(2.1) 
(2.2) 




(2.3) 
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Notice that condition (HI) implies the strong Legendre condition along an 

extremal pair {u{t),Xt): 

{Xt){v,v) < -a\v\'^, t e [0,ti], V e M™, a > 0, 

u—u{t) 

i.e., the extremal At is regular [2]. 

Moreover, we assume also that the following condition is satisfied: 

(H2) For any A G T* M , the function u /i„(A), u G U , has a maximum point 
u{X) e U: 

hu(x)W = max/i„(A), A e T*M. 

In terms of work [1], condition (H2) means that T*M is a regular domain of the 

Hamiltonian ^.^(A). Condition (HI) means that the function u i— !■ hu{X) has no 
maximum points in addition to w(A). At the maximum point l^^^-ji^-) (A) = 
for all A S T*M. By implicit function theorem, the mapping A i-^ 7i(A) is 
analytic. The maximized Hamiltonian H(X) = /Iu(a)(A), A € T* M . is also 
analytic. The extremal Xt is a trajectory of the corresponding Hamiltonian 
vector field: At = H{Xt), and the extremal control is u{t) = u{Xt). 

2.2 Second veiriation and its Morse index 

Consider the endpoint mapping for problem (2.1)-(2.3): 

Ft:U = L^{[0,t],U) ^ M, n(-) ^ (g„(t), J*M), (2.4) 

where g„(-) is the trajectory of the control system (2.1) with the initial condition 

9ii(0) = (Jo corresponding to the control u = u(-). Since m G W is an extremal 
control, it follows that the differential {first variation) : T^U — > Tg^^^^M is 
degenerate, i.e., not surjective, for all t G {0,ti], see [2]. 
Introduce one more important hypothesis: 

(H3) The extremal control u{-) is a corank one critical point of the endpoint 
mapping Ft, i.e., 

codim Im i:)5Ft = 1, t& (0,ti]. 
Condition (H3) means that there exists a unique, up to a nonzero factor, ex- 
tremal At corresponding to the extremal control u{t). 

For any extremal control u eU there is a well-defined Hessian {second vari- 
ation) [2] of the endpoint mapping a quadratic mapping 

Hess„ft : KerD„Ft ^ Coker L»„i^t = Tg^(t)M/ Im D^Ft. 

Condition (H3) means that dim (r,^^(t)M/ImDsFt) = 1 for all t G {0,ti], thus 
uniquely, up to a positive factor, there is defined a quadratic form 

Qt = At Hesse Ft : Ker D^Ft ^ K, t e {0,ti], (2.5) 
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the projection of the second variation to the extremal At. 

The Morse index of a quadratic form Q defined in a Banach space £ is the 
maximal dimension of the negative space of the form Q: 

indQ = maxjdimL \ L C JC, Q\l\^o} < ^l- 

The kernel of a quadratic form Q{x) is the space 

Ker Q = {x e £ I Q{x, y) =0\/ y € £}, 

where Q{x, y) is the symmetric bilinear form corresponding to the quadratic 
form Q{x). A quadratic form is called degenerate if it has a nonzero kernel. 
The multiplicity of degeneration of a form Q is equal to dimension of its kernel: 
dgn Q = dim Ker Q. 

Now we return to the quadratic form Qt given by (2.5) — the second varia- 
tion of the endpoint mapping for the extremal pair u{t) , At of the optimal control 
problem (2.1)-(2.3). We continue the quadratic form Qt from the space Loo to 
the space L2 by continuity, and denote by Kt the closure of the space Ker D^Ft 
in L2[0,t]. 

Proposition 2.1 (Propos. 20.2 [2], Th. 1 [12]). Under hypotheses (HI) , 
(H3) , the quadratic form Qtlxt positive for small t > 0. In particular, 
ind Qt\xt = Z^*" small t > 0. 

An instant € (0, ti] is called a conjugate time (for the initial instant t = 0) 
along the extremal At if the quadratic form Qt. is degenerate. In this case 
the point qu{t*) = 7r(At,) is called conjugate for the initial point qo along the 
extremal trajectory qu{-)- 

Proposition 2.2 (Th. 1 [12]). Under hypotheses (HI) , (H3) : 

(1) conjugate points along the extremal At are isolated: Q <t\ < ■■ ■ < t^ < 

h, 

(2) Morse index of the second variation is expressed by the formula 

ind gtlx, = ^{dgnQti | < < i}. 

Local optimality of extremal trajectories is characterized in terms of conju- 
gate points. Speaking about local optimality of extremal trajectories in calculus 
of variations and optimal control, one distinguishes strong optimality (in the 
norm of the space C([0, ti],M)), and weak optimality (in the norm of the space 
C^{[0,ti], M)). Under hypotheses (HI) -(H3) , normal extremal trajectories 
lose their local optimality (both strong and weak) at the first conjugate point [2]. 
Thus in the sequel, when speaking about local optimality, we mean both strong 
and weak optimality. 

Proposition 2.3 (Propos. 21.2, Th. 21.3 [2]). Let conditions (HI) -(H3) be 
satisfied. 
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(1) If the interval (0,ti] does not contain conjugate points, then the extremal 
trajectory q{t), t G [0,ti], is locally optimal. 

(2) If the interval (0,ii) contains a conjugate point, then the extremal trajec- 
tory q{t), t € [0,fi], is not locally optimal. 

2.3 Exponential mapping 

We will add to hypotheses (HI) -(H3) one more condition: 

(H4) All trajectories of the Hamiltonian vector field H{X), A G T*M, are con- 
tinued to the segment t G [0./:i]. 

Consider the exponential mapping for time t: 

Exp, : N = T;^M ^ M, Exp,(A) = tt o e*^(A) = q{t), t G [0, ti]. 

One can construct a theory of conjugate points in terms of the family of the 
subspaces 

A{t) = e-'^T^^{T;^,^M)CT^„{N), 

via linearization of the flow of the Hamiltonian vector fleld H along the ex- 
tremal At. 

2.4 Maslov index of a curve in Lagrange Grassmanian 

First we recall some basic facts of symplectic geometry, see details in works [1,4]. 
Let I],o- be a symplectic space, i.e., S is a 2n-dimensional linear space, and a 
is a nondegenerate skew-symmetric bilinear form on E. The skew-orthogonal 
complement to a subspace F c S is the subspace T"^ — {x G T, \ a{x, F) = 0}. 
Since a is nondegenerate, it follows that dimF^ = 2n — dimF. A subspace 
F c E is called Lagrangian if F = F^, in this case dimF = n. The set of 
all Lagrangian subspaces in S is called Lagrange Grassmanian and is denoted 
as L{Ti), it is a smooth manifold of dimension n(n + l)/2 in the Grassmanian 
G„(E) of all n-dimensional subspaces in E. 

Fix an clement H e Define an open set = {A e L(S) | Ann = 0}. 

The subset Mn = L{T.) \ = {A e L(E) | A n H ^ 0} is called the train 
for n. The set Mn is not a smooth submanifold in LiTi), but it is represented 
as a union of smooth strata: Mn = Ufc>i M^^\ where m'^^^ = {A e L(S) | 
dim(A n n) = A:} is a smooth submanifold of i(S) of codimension k(k + l)/2. 

Consider a smooth curve A(<) e t G [to, ii], i-e., a family of Lagrangian 

subspaces in E smoothly depending on t. Suppose that A(io),A(ti) G H*. 
Maslov index /xn(A(-)) of the curve A(-) is the intersection index of this curve 
with the set A^n- 

In greater detail, let the curve A(-) do not intersect with Mn \ -'^n'', this 
can always be achieved by a small perturbation of this curve. For the smooth 
hypersurface A^n"* C L{Yi), one can define its coorientation in an invariant way 
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as follows. Any tangent vector to at a point A € can naturally be 

identified with a certain quadratic form on A. Take a tangent vector A(t) G 
Ta(()L(E) to a smooth curve A{t) G LCE). Choose a point x G A{t) of the n- 
dimensional space A(i) C E. Choose any smooth curve r x(t) in E such that 
x{t) G A(t) for all r, and x{t) = x. Then the quadratic form A{t){x), x G A{t), 
is defined by the formula A{t){x) = a{x,x{t)). One can show that a{x,x{t)) 
docs not depend upon the choice of the curve x(t), i.e., one obtains a well- 
defined quadratic form A{t) on the space A{t). Moreover, the correspondence 
A I— > A, A e TaZ/(E), defines an isomorphism of the tangent space TaL{T,) and 
the linear space of quadratic forms on A, sec [1]. 

Maslov index yL(n(A(-)) is defined as the number of transitions of the curve 
A(-) from the negative side of the manifold (i.e., with A{t) > 0) minus 

number of reverse transitions (with A{t) < 0), taking into account multiplicity. 

The fundamental property of Maslov index is its homotopy invariance [3]: for 
any homotopy A%t), t € s G [0,1], such that A" {I^q) , A^tf) G for all 

s G [0, 1], we have /Un(A"(-)) = ^n(A^(-)). This fact is proved in the same way 
as homotopy invariance of the usual intersection index of a curve with smooth 
cooriented surface. 

For monotone curves in Lagrange Grassmanian i(E) there is the following 

way of evaluation of Maslov index. 

Proposition 2.4 (Cor. I.l [1]). Let A{t) <0, t e [ta,ti], and let {t G [io,ii] | 
A(t) n n 7^ 0} be a finite subset of the open interval (fo, ti). Then 

Mn(A(-)) = - J2 dim(A(t)nn). (2.6) 

In fact, in Cor. I.l [1], there is given a statement for a nondecreasing curve 
{A{t) > 0), then in the right-hand side of formula (2.6) the sign minus is absent. 
As indicated in the remark after Cor. I.l [1], the passage from nondecreasing 
curves to nonincreasing ones is obtained by the inversion of direction of time 
t to + ti — t. 

The theory of Maslov index can be used for computation of Morse index for 
regular extremals in optimal control problems. 

2.5 Morse index and Mcislov index 

let At, t G [0,/:i], be a normal extremal of the optimal control problem (2.1)- 
(2.3), and let the hypotheses (HI) -(H4) be satisfied. Consider the family of 
quadratic forms Qt given by (2.5). 

The extremal At determines a smooth curve 

A{t) = e-'"Tx,{T;(,^M) G i(E), t G [0,ti], 

in the Lagrange Grassmanian i(E), where E = Txg{T* M). The initial point of 
this curve is the tangent space to the fiber A(0) = 11 = T\g{T*^M). The strong 
Legendre condition (see (HI) ) implies monotone decreasing of the curve A{t): 



8 



the quadratic forms A(t) < 0, t G [0,ti], see Lemma 1.4 [1], thus its Maslov 
index can be computed via Propos. 2.4. 

On the other hand, the following important statement establishes relation be- 
tween Morse index of the second variation Qt and Maslov index of the curve A(t). 

Proposition 2.5 (Th. 1.3, Cor. 1.2 [1]). Let hypotheses (HI) -(H4) be satisfied. 
Then: 

(1) An instant t G {0,ti\ is a conjugate time iff A{t) fl 11 ^ 0. 

(2) If A{ti) n n = 0, then there exists i > such that 

ind Qt,\j,^^ = -/zn(A(-)l[to,ti]) ^ to e (0,i). 

(3) // {t G (0, ti] I A{t) n n ^ 0} is a finite subset of the open interval (0, ti), 
then 

indQtJ^^^= J2 dim(A(t)nA(0)). 

fe(o,ti) 

Item (1) of Propos. 2.5 implies obviously the following statement. 

Corollary 2.1. Let hypotheses (HI) (H4) hold. An instant t € (0,<:i) is a 

conjugate time iff the mapping Exp^ is degenerate. 

Proof The condition A{t) nU^O means that e*^(n) n Tx, {T*^t)^) ^ 0, which 
is equivalent to degeneracy of the mapping Exp^ = tt o e*^. □ 

Due to Propos. 2.5, wc obtain a statement on homotopy invariance of Maslov 

index of the second variation. 

Proposition 2.6. Let (u*(f),A|), t G [0,ff], s G [0,1], be a continuous in pa- 
rameter s family of normal extremal pairs in the optimal control problem (2.1)- 
(2.3) satisfying the conditions (HI) -(H4) . Assume that for any s G [0, 1] the 
terminal instant t = tl is not a conjugate time along the extremal . Then 



ind Qf 



= ind QfO 



(2.7) 



Proof. It follows from continuity and strict monotonicity of the curves A* (t) = 
e^*^TAj(r*,(j)M), q^l^t) = 7r(Af) that there exists i > such that t < ts for all 
s G [0, 1] and any instant t E (0, t) is not a conjugate time along the extremal Xf. 
According to item (2) of Propos. 2.5, we have 

ii^dQtflxts =-MA^(-)l[t„t3]), Vtoe(0,i) VsG[0,l]. (2.8) 

For all s G [0, 1] we have A*(fo)nn = A*(tf)nn = 0. Then homotopy invariance 
of Maslov index imphes that the function s i-^ A*n( A*(-)|j(^ j^j) is constant at the 
segment s G [0,1]. Thus equality (2.8) implies the required equality (2.7). □ 
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The following statements can be useful for the proof of absence of conjugate 

points by homotopy or limit passage. 

Corollary 2.2. Let all hypotheses of Propos. 2.6 be satisfied. If an extremal 
trajectory — 7r(At), t G {0,t1], does not contain conjugate points, then the 
extremal trajectory q^{t) = 7r(Aj), t G (0,f}], does not contain conjugate points 
as well. 

Proof. A regular extremal does not contain conjugate points iff its Maslov index 
is zero, so the statement follows from Propos. 2.6. □ 

Corollary 2.3. Let (u*(f),Af), t € [0,+oo), s G [0,1], be a continuous in pa- 
rameter s family of normal extremal pairs in the optimal control problem (2.1)- 
(2.3) satisfying the hypotheses (HI) -(H4) . Let for any s G [0, 1] and T > 
the extremal A| have no conjugate points for t G (0, T]. Then for any T > the 
extremal Aj has no conjugate points for t G (0, T] as well. 

Proof. Fix any T > 0. By Propos. 2.2, conjugate points along the extremal Aj 

are isolated, thus there exists an instant ti > T that is not a conjugate time 
along Aj. Consider the family of extremals A|, t G [0,ti], s G [0, 1]. Corollary 2.2 
implies that the extremal A^ has no conjugate points for t G {0,ti], thus for 
t G (0, T] as well. □ 

2.6 Preliminary remarks on Euler's problem 

In this subsection we show that Euler's elastic problem satisfies all hypotheses 
required for the general theory of conjugate points described in Subsecs. 2.1-2.5. 
Recall [11] that Euler's problem is stated as follows: 

q = Xi{q) +uX2{q), qe M = R'^ x S\ u G K, (2.9) 
q{0) = 90, q{ti) = qi, ii fixed, (2.10) 

1 Z'*^ 

J = - / u^dt min, (2.11) 

2 Jo 

where = cos + sin , ^2 - [Xi, X2] = = sin - cos J^. 

This problem has the form (2.1)-(2.3), and the regularity conditions for M, 
/, (p are satisfied. 

In terms of the Hamiltonians h,{X) = {\,X,), A G T*M, i = 1,2,3, the 
normal Hamiltonian of PMP for Euler's problem is hu{X) = hi{X)+uh2{X) — ^u^. 
d^hu 

We have ^ " = — 1 < 0, i.e., hypothesis (HI) holds. 

a V, 

Condition (H2) obviously holds. 

Let u{t) be a normal extremal control in Euler's problem. Corank of the 
control u{t) is equal to dimension of the space of solutions to the linear Hamil- 
tonian system of PMP Aj = hi{Xt) + u{t)h2{Xt), i.e., to the number of distinct 
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nonzero solutions to the Hamiltonian system corresponding to the maximized 
Hamiltonian H = hx + }^h%: 

At = hi{\t) + h2h2{Xt), u{t) = h2{Xt). (2.12) 

We are interested in the number of distinct nonzero solutions to the vertical 
subsystem of system (2.12): 

hi = -h2h3, f/3=c, 

/i2 = /13, ^ <c=-rsin/J, (2.13) 

hz = hihz [»^ = 0, 

where hi = — rcos/5, /12 = c, ~ — rsin/3, sec [11]. 

To the extremal control u{t) = 0, there correspond two distinct nonzero 
extremals (/ii, /i2, ^3)(At) = (±r, 0,0), r ^ 0, so in this case coranku = 2. 

If u(t) ^ 0, then Ci = u{t) ^ 0. Then the function Ci determines uniquely 
via system (2.13) the functions rsin/Jj = —ct and rcos/Jj = —Ct/ct- So the 
curve {hi,h2,hz){Xt) ^ is uniquely determined. Consequently, coranku = 1 
in the case u{t) ^ 0. 

Notice that the control u{t) = is optimal, thus in the sequel in the study 
of optimality of extremal controls we can assume that their corank is equal to 1, 
i.e., hypothesis (H3) is satisfied. 

Finally, hypothesis (H4) is also satisfied since the Hamiltonian field H is 
complete (its trajectories are parametrized by Jacobi's functions determined for 
all t e M). 

Summing up, all hypotheses (HI) -(H4) are satisfied for Euler's elastic 
problem, so the theory of conjugate points stated in this section is applicable. 



3 Conjugate points on inflectional elasticae 

In this section we describe conjugate points on inflectional elasticae in Eu- 
ler's problem. We perform explicit computations and estimates on the basis 
of parametrization of extremal trajectories obtained in [11]. 

We base upon the decomposition of the preimage of the exponential map- 
ping T*^M = N = uJ^iNi introduced in [11]. In this section we consider the 
case A € Ni. In Subsec. 8.2 [11] was obtained a parametrization of the expo- 
nential mapping in Euler's problem Exp^ : {(p,k,r) 1-^ {xt,yt,Ot) in terms of 
elliptic coordinates in the domain Ni. By virtue of Corollary 2.1, an instant 
f is a conjugate time iff the mapping Expj is degenerate, i.e., iff its Jacobian 

J = vanishes. A direct computation using parametrization of the 

d{ip,k,r) 



11 



exponential mapping obtained in Subsec. 8.2 [11], yields the following: 



d{xt,yt,et) ^ 1 d{xt,yusin{9t/2)) ^ 32fc 

d{ip,k,r) V^cos(6't/2) d{ip,k,^) (1 - fc2)r3/2A2 ^' 

(3.1) 

= ao + aiz + 02-2^, ^; = sn^r e [0, 1], (3.2) 

a'2 = —k^ snp xi, (3.3) 

012 + ai + ao = (1 - fc^) snp a;i, (3.4) 

ao = fi{p,k)x2, (3.5) 



a;i = — dnjj(2snp dnp E'^(p) + ((4fc2 — 5)psnp dnp 

+ cnp(3-6fc2sn2p))E2(p) + {{ik^ - 5) cup {1 - 2k^ sri^ p) p 
+ snp dnp(4p^ - 1 + fc^(6sn^p - 4 - 4p^)))E(p) 
+ p snp dnp(l - (1 - fc2)p2 + fc2(^4^2 _ 5)sn2^^ 

+ 2cnp(/c2sn2p dn^p + (1 - fc^)(l - "^k^ sn^ p)p^))), (3.6) 
= cnp(2(l-fc2)E(p) -E2(p) - (I - p)p^) 

+ snp dnp (E(p) - (1 - k^)p), (3.7) 
fiip,k) = snp dnp - (2E(p) -p)cnp, (3.8) 
p=^t/2, r = \/r(v + t/2), A = 1 - fc^ sn^p sn^ r . 

Here cn, sn, dn, E are Jacobi's functions, see details in [11]. 
3.1 Preliminciry lemmcis 

In this subsection wc describe roots and signs of the functions ciq and aa + ai + ap 
that essentially evaluate the numerator of Jacobian J at the extreme points 
z = and 1 respectively, see (3.1), (3.2). 

3.1.1 Roots of the function qq 

Roots of the function /i(p) defined in (3.8) where described in work [10]. For 
completeness, we cite the statements we will need in the sequel. 

Proposition 3.1 (Lemma 2.1 [10]). The equation 2E{k)-K{k) = 0, A; G [0, 1), 
has a unique root ko € (0, 1). Moreover, 

k e [0, ko) 2E-K> 0, 

fce(fco,l) 2E-K<Q. 

Here and below K{k) and E{k) are complete elliptic integrals of the first 
and second kinds, see [8,11,13]. Numerical computations give the approximate 
value fco « 0.909. 
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Proposition 3.2 (Propos. 2.1 [10]). For any k G [0,1) the function fi{p,k) 
has a countable number of roots p\, n e Z, localized as follows: Pq = ^ and 

pI e {-K + 2Kn, K + 2Kn), n e Z. 

Moreover, for n e N 

fc e [0, fco) ^ pi€ {2Kn, K + 2Kn) , 

k = ko => Pn = 2Kn, 

fce(fco,l) ^ p\€{-K + 2Kn,2Kn), 

where ko is the unique root of the equation 2E{k) — K{k) = 0, see Propos. 3.1. 

Now we establish the signs of the function fi{p) between its zeros pj^. 
Lemma 3.1. For any m = 0,1,2, ... , we have: 

P e {P2 P2m+1 

P e {plm+l,plm+2) Mp) < 0. 

Proof. By virtue of the equality 

/i (p) \ ' sv? p dv? p 



cnp J cn^p ' 

the function /i {p) /cup increases at the segments of the form [—K + 2Kn, K + 
2Kn], n e Z. So the function /i(p)/cnp, as well as fi{p) changes its sign at the 
points e {—K + 2Kn, K + 2Kn). It remains to verify that fi{p) is positive 
at the first interval {pl,p\) = (0,p1). We have /i(p) = //3 + o(p^) > 0, p ^ 0, 
and the statement follows. □ 

Now we describe zeros of the function X2 that enters factorization (3.5) of 
the function ao- 

Lemma 3.2. The function X2{p) given by (3.7) has a countable number of roots 
P = Pn^ > 0. We have p^^ = and p^^ € {2Kn, K + 2Kn) for n e N, moreover, 

k<ko p^' e{pl,K + 2Kn). (3.9) 

Further, 

P&{plLA+i) ^2{P)>0, (3.10) 

P&{PTm+l>PTm+2) X2{p)<0, m = 0,1,2, ... . (3.11) 

Proof. First we show that the function 

X2{P) 



snp dnp 



increases when p G {2Kn, 2K + 2Kn). (3.12) 



13 



Direct computation gives 

^)'=^^, (3.13) 
snp dnp J sn^p dn p 

X3 = e{cn^p E(p) + + (1 - k^){E{p) + pf, (3.14) 
Q = (1 — k^)psv? p — cap snp dnp, (3 = —pdv^p. 

Since 

Hp) + /3 = ^A:V + o(p'), cn^p E(p) + a = ^(1 - &2)p3 + ^(^3)^ 

then E(]3) + /3 ^ 0, cn^p E(p) + a ^ 0. So the function x^{p) given by (3.14) is 
nonnegative and vanishes only at isolated points. By virtue of equality (3.13), 
assertion (3.12) follows. 
Further, we have 

^■2\p=2Kn = cnpa;4(p), 

X4, = -((1 - fc^)(E(j5) - pf + fc^ e2(p)) < for all p ^ 0. 

Thus 

p=2K + AKn => cnp < 0, X2 > 0, 
p = AKn => cnp > 0, X2 < 0. 

Consequently. .T2/(snp dnp) ±oo as p — > 2ii'n =F 0, n € N. Moreover, it 
follows from the asymptotics 

X2 = (4/45)fe2(l-fc2)/ + o(/), p^O, (3.15) 

that a;2/(snp dnp) ^ +0 as p +0. 
Thus 

pe(0,2K) ^ ^5 >0, 

snp dnp 

p e (2iirn, 2iir + 2Kn) => increases from — oo to + oo. 

snp dnp 

So there exists a unique root of X2(p)/(snp dnp), thus of a;2(p) at the interval 
{2Kn, 2K + 2Kn), we denote it as p^\ 

Now we localize p^^ w.r.t. the point K + 2Kn. We have 

x^\p=K+2Kn = dnp snp(E(p) - (1 - k^)p), 

E(p) - (1 - k^)p = fc^ / cn^tdt>0, p> 0. 
Jo 



Now 



p = K + AKn snp = 1, a;2 > > 0, 

snp 

p = 2K + AKn snp = -1, 0:2 < > 0. 

snp 
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Consequently, pf^^ g (^2Kn, K + 2Kn) for all n G N. 

Let k < ko, thenp^ G {2Kn, K+2Kn), we clarify now the mutual disposition 
of the points and p^^ in this case. By virtue of (3.8), 

/i(p) = E(p) = (dnp snp/cnp +p)/2. 

Direct computation gives 

Ie(p)=( d„ p s„ p / c„ p +p)/2 = 8 cn^ p ( sn p dn p - p cn p ) E(p) . 

Since for p = pjj we have cnp ^0, it follows that for p = p^ the functions 
X2 and snp dnp — pcnp have the same sign. Then both for p = pli_i G 
{4KI - 2K, AKl - K) and for p = p^; G {AKl, AKl + K) wc obtain X2/ snp < 0. 
Consequently, p\ < p^ for all n G N, i.e., inclusion (3.9) is proved. The roots 
p^ are localized as required. 

For p > the functions X2 and snp have distinct roots, so it follows 
from (3.12) that X2 changes its sign at the points p^ , n G N. The distri- 
bution of signs (3.10), (3.11) follows from the fact that at the first interval 
(Po^jPD = (OjPD tlie function X2 is positive, see (3.15). □ 

For p > 0, the function ao vanishes at the points p = p^ and p = p^ defined 
and localized in Propos. 3.2 and Lemma 3.2. Now decomposition (3.5) and 
Lemmas 3.1, 3.2 imply the following statement on distribution of signs of the 
function ao- 

Lemma 3.3. Let k G (0,1). If p G (0,p}), then oq > 0. For any n G N, if 
p G (p^,P^^), then ao < 0, and if p G (p^^pJ^_,_l), then oq > 0. 

3.1.2 Roots of the function oq + ai + 02 

In order to obtain a similar description for the function 00 + 01 + 02, we have 
to describe roots of the function Xi, see decomposition (3.4). 

Lemma 3.4. For p > 0, the function xi{p) defined by (3.6) has a countable 
number of roots po = 0, pf^^ G {Pn,Pn+i), n gN. Moreover, 

P^ivVraA+i) ^ a;i(p)>0, (3.16) 

'P2m+2) ^ ^1 

(p) < 0. (3.17) 

Proof. Direct computation gives 

( xx{p) y_ xr,{p) , . 

\Anpf,{p)) - AflipY ^'-'"^ 
xr, = k^{cnpE4p + af + (1 - k'^){pE2+(3f > 0, 

E2 = 2E(p) -p, E4 = cnp(2E(p) -p) - 2 snp dnp, 
a = (1 + sn^ p - 2fc^ sn^ p )£;| + 4 cnp snp dnp (1 - 2fc^)£;2 

+ 4(2A:2 - l)sn2p dn^p, 
/? = (2fc^sn^p - l)i;| +8fc^cnp snp dnp £2 - Sfc^sn^p dn^p. 
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Since 

45 

the function 2:5 (p) is nonncgativc and vanishes at isolated points. In view of 
equality (3.18), the function Xi{p) / {dnp fi{p)) increases at the intervals where 

Now we find the sign of xi at the points pj^. We have 
I _ xejp) 

^llE(p) = (dnp snp/cnp+p)/2 4 cn^ p ' 

Xq{p) = x^ + xIp + xIp"^, 

Xq{p) = — cn^ p dnp (1 — sn^ P (2 — sn^ p )), 

xl{p) = 2cnp snp (1 — 2fc® sn® p — sn'^p (3 + sn^p) + sn'^p (4 + sn^p)), 
Xg(p) = — dn^p sn^p(l — A:^sn^p(2 — sn^p)), 

notice that 1 — A;^sn^p(2 — sn^p) = dn^p + A;^(l — fc^)sn'^p > 0. Consider 
the discriminant 

X6d = {xlf -4x^x1 = -16A:2(1 - A;^) cn^p sn*^ jjdn^p 

of the quadratic polynomial xc,{p). If A; 7^ fco, then for p = p\we have cnp ^0, 
snp ^ 0, thus x1,x1,XQd < and < 0. \i k = fco, then for p = p\ we have 
cnp ^ 0, snp 7^ 0, thus x% = 0, x| < 0, x^d = 0, = and xq = x^p^ < 0. 

So for all k G (0,1) if p = p\ > 0, then sgnxi = — sgncnp. If p = 
P21-1 £ i'^K^ - '2K,AKl - K), then cnp < 0, thus xi > 0. Similarly, if 
p = pij G {4:Kl, 4KI + K), then cnp > 0, thus Xi < 0. 

Consequently, 

p G (0,p}) =J> '''^^fl increases from to +00 => xUp) > 0, (3.19) 
dnp/i(p) 

p G (pi,pi_i_i), n G N '''"^^f-^ increases from — 00 to +00, 

Kfn'fn+Lj' dnp/i(p) 

thus xi has a unique root p^^ G (p^,p^j_|_i). 

The required signs of the function .Ti(p) at the intervals (3.16), (3.17) fol- 
low from the inequality at the first interval (3.19), and from the fact that 
Xi(p)/(dnp/i(p)), thus Xi{p) changes its sign at the points p^\ n G N. □ 

Remark. By virtue of decomposition (3.4), we have the equality 

{p > I ao + ai + 02 = 0} = {p > I snp = 0} U {p > I xi = 0} 

= {2Km I m G N} U {p^' | n G N}. (3.20) 
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In order to have a complete description of roots of the function 00 + 01 + 02, 
one should describe mutual disposition of the points 2Km and . Numerical 
computations show that some of these points may coincide one with another. 
For example, numerical computations yield the following relations between the 
first roots in the families (3.20): if /c G (0, A;), then p^^ > 2K; if k = k, then 
= 2K; if e [k, 1), then < 2K for a number k w 0.998. We do not go 
into details of this analysis, but in the sequel we allow different possibilities of 
mutual disposition of the roots and 2Km. 

3.2 Bounds of conjugate time 

In this subsection we estimate the first conjugate time in Euler's problem along 
inflectional elasticae. 

We obtain from equalities (3.3), (3.4) that 02 = — fc^/(l — A:^)(ao + oi + 02), 
thus the Jacobian appearing in (3.1), (3.2) can be represented as 



notice that (1 — k'^z) /(I — k"^) > 0. In order to describe the first conjugate point 
along an extremal trajectory q{t) = w o e*^(A), X G Ni, it suffices to describe 
the first positive root of the function Ji for fixed k, z: 



This minimum exists since by virtue of regularity of normal extremals, small 
intervals p e (0,e) do not contain conjugate points. Below in the proof of 
Th. 3.1 we show this independently on the basis of explicit expression for the 
function Ji. 

Theorem 3.1. Let X € Ni. For all k e (0,1) and z G [0,1], the number 
p'l°'^\k,z) belongs to the segment bounded by the points 2K{k), Pi{k), namely: 



(1) fce(0,fco) pr^€[2K,pl], 

(2) k = ko = = Pi, 

(3) fce(fco,i) pr^e\pl,2K]. 

Moreover, for any k G (0, 1) there exists s = e{k) > such that: 
{!') IfkG (0,fco), then 



Ji{p,k,z) = (1 - z)ao + z{l - k'^z)/{l - k^){ao + ai + 02) 



(3.21) 



p^°°J(fc, z) = min{p > | Ji(p, k, z) = 0}. 



p G (0, 2K) ^ Ji > 0, 



(3.22) 
(3.23) 



(2') Ifk = ko, then 



pe{0,2K) ^ Ji>0, 

p G {2K, 2K + e) ^ Ji < 0, 



(3.24) 
(3.25) 
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(3') Ifke (A;o,l), then 

p€{0,pl) Ji>0, (3.26) 

moreover, 

(3' a) in the case p'^^ € {Pi,'2K) : 

P^iPi\p?+e) Ji<0, (3.27) 

(3'6) in the case = 2K : 

p = 2K ^ ^ Ji < 0, (3.28) 

(3'c) in the case p^' e (2-ftr,pf) : 

p e (2ii', 2iV' + £) Ji < 0. (3.29) 

Proof. It is easy to see that by virtue of continuity of the function Ji(p), items 
(1'), (2'), (3') imply respectively items (1), (2), (3), so we prove statements (1'), 
(2'), (3'). 

(1') Fix any k G (0, ko), then 2K < pi, see Propos. 3.2. 
If p e {0,2K), then Lemmas 3.1, 3.2, 3.4 and decompositions (3.5), (3.4) 
imply the following: 

/i > and X2 > =^ oo > 0, 

snp > and a;i > => ao + oi + 02 > 0. 

Then representation (3.21) gives the inequality Ji(p, 2;) > for all 2; € [0, 1] and 
allpG {0,2K). Implication (3.22) follows. 

Lemmas 3.2, 3.4 imply respectively that Pj^ e {p\,3K), pj^ e (Pi,P2)- 
Denote pi = min(p'^^,p^^) > p\. 

If p G (Pi,Pi), then we obtain from Lemmas 3.1, 3.2, 3.4 and decomposi- 
tions (3.5), (3.4) the following: 

/i > and X2 > ^ ao < 0, 

snp < and a;i > ao + a\ + < 0. 

Representation (3.21) implies that Ji(p, z) < for all z G [0,1] and all p G 
(p},pi), i.e., implication (3.23) is proved for s = pi — pi > 0. 
(2') Let k = ko. Similarly to item (1'), 

p G (0, 2K) ao > and oq + ai + 02 > => Ji > 0, 

pG(2ii',pi) ao < and ao + ai + a2 < Ji < 0, 

where pi = min(p5^i ,p^=) > 2K. Thus implications (3.24), (3.25) follow for 

e = Pi - > 0. 

(3') Let fc G (/Co, 1), then pl{k) < 2K{k). 
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Let p e {0,Pi)- Then we have the following: 

/i > and X2> ao > 0, 

snp > and a;i > => ao + ai + 02 = 0. 

Thus Ji > 0, and implication (3.26) is proved. 

(3'a) Consider the case p^' € {pl,2K). Let p G (p^' , 2K), then: since /i < 
and X2 > 0, then ao < 0; since snp > and Xi < 0, then ao + ai + 02 < 0. 
Thus Ji < 0, and implication (3.27) follows for e = 2K — p^^ > 0. In this case 

pr\z) e bl.Pr] c \pl,2K) yze[o, 1]. (3.30) 

(3'6) Consider the case p^^ = 2K. Let p = 2K, then: since /i < and 
X2 > 0, then ao < 0; since sup — xi = 0, then ao + ai + 02 = 0. Consequently, 
Ji < 0, and implication (3.28) follows. 

(3'c) Finally, consider the case p^' e {2K,pl). Let p e {2K,Tain{pl\pl-')), 
then: since /i < and X2 > 0, then qq < 0; since snp < and xi > 0, 
then ao + ai + a2 < 0. Thus Ji < 0, and implication (3.29) is proved for 
e = min(pf ,pf)-2J!r>0. □ 

Remark. As one can see from inclusion (3.30), for p^^ e {pl,2K) the range of 
the function p'{°'^'{z) is strictly less than the segment [p|,2i>l']. Judging by the 
plots of the function p1°'^\z), z = sn^ t , this function is smooth and strictly 
monotone at the segment r G [0, JsT], see Figs. 1-4. 




Figure 1: p = p'^'^\k,T), k G {0,ko) Figure 2: p = (fc, r) k = ko 

From decompositions (3.5), (3.4) and Lemmas 3.2, 3.4 we obtain the fol- 
lowing description of all (not only the first) conjugate points for the cases 
z = sn^T = or 1 (i.e., for elasticae centered respectively at its vertex or 
inflection point). 

Corollary 3.1. Let X € Ni and k e (0, 1). 

(1) Ifz = 0, then {p>0\ Ji{p,z) = 0} = K | n e N} U ^ | m e N}. 

(2) Ifz = l, then {p>0\ Ji{p,z) = 0} = {2Kn \ nen}U{p^ \ m€ N}. 
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2K 




K 



Figure 3: p = p^°^^{k,T), k G 
{ko,l),2K<pl^ 



Figure 4: p = pl°'^\k,T), k G 
(fco.l), 2K>pl^ 



Remark. According to Lemma 3.2 and Propos. 3.2, in item (1) of Cor. 3.1 all 
the roots pj, and p^ are pairwise distinct. Although, in item (2) some of the 
roots 2Kn and p^ may coincide one with another, see the remark at the end 
of Subsec. 3.1. 

Now we apply preceding results in order to bound the first conjugate time 
along normal extremal trajectories in the case X G Ni: 

t'^'^\X) = mm{t > \t conjugate time along trajectory q{s) = Exp^(A)}. 

Theorem 3.2. Let A = (fc, ip, r) G Ni. Then the number t'{°^\X) belongs to the 

. , , , . 4K(k) 2p]{k) 
segment with the enapomts — ^— , — ^— , namely: 

\/r Jr 



(1) fee (o,fco) 

(2) /e = A:o = 

(3) k G (fco, 1) 



^conj g 



4K{k) 2p\(k) 



j _ 4K{k) _ 2p\(k) 

II — ^ — r- 



r \/r 

2p\{k) AK{k) 



Proof. By Corollary 2.1, an instant i > is a conjugate time iff 
d{xt,yt,Ot) 32k 



J{t,k,ip,r) = 



d{(p,k,r) (1 - A:2)r3/2A2 
p=V^t/2, T = y/r{ip + t/2), 2;=snV 

see (3.1). 



Ji{p,k,z) = 0, 

A = 1 — A;^ sn^ p sn^ r , 



(1) Let k G {0,ko), then i^M < ^^iM. According to item (1') of Th. 3.1, 



\/r -v/r 
for some e = s{k) > we obtain the chains: 



t e 0, 



7r 



PG{0,2K) 



J{t, k,(p,r) < \/ ip, r, 
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i'2p\(k) 2(p\(k) + e)\ , ,,,, 



/r wr 



J{t,k,(p,r) > V iy5,r. 



By virtue of continuity of the function J w.r.t. t, we obtain the required inclusion 



^conj g 



AK{k) 2pl{k) 



/r 

Statements (2), (3) of this theorem follow similarly from items (2'), (3') of 
Th. 3.1. □ 

In Section 12 of work [11], was defined a function t : N ^ {0,+oo] that 
provides an upper bound on cut time in Euler's elastic problem, see Th. 12.1 [11]. 
It follows from formula (12.2) [11] that 

Comparing this equality with Th. 3.2, we obtain the following statement. 

Corollary 3.2. IfXeNi, then fi°"j(A) > t(A). 

A natural measure of time along extremal trajectories in Euler's problem is 
the period of the pendulum T{k) = 4:K{k) j \fr. In terms of this measure, the 
bounds of Th. 3.2 are rewritten as follows. 

Corollary 3.3. Let Xg Ni. Then: 

(1) k e (0, ko) => tr^ e [T, tl] C [T, 3T/2), tl = 2pl/^ e (T, 3r/2), 

(2) k = ko ^ tr'' = T, 

(3) k G {ko, 1) ^ tr'' e [tl,T] C (T/2,T], tl = 2p\l^r e {T/2,T). 

It is instructive to state the conditions of local optimality for elastica in 
terms of their inflection points. 

Corollary 3.4. Let A € A''!, and let T = {7^ = {xs,ys) \ s e [0,t]}, q{s) = 
{xs,ys,Os) = Exp(As), be the corresponding elastica. 

(1) // the arc T does not contain inflection points, then it is locally optimal. 

(2) If k e (0, ko] and the arc T contains exactly one inflection point, then it 
is locally optimal. 

(3) // the arc T contains not less than three inflection points in its interior, 
then it is not locally optimal. 

Proof. (1) If the elastic arc T does not contain inflection points, then its cur- 
vature Cg = 2k^/rcn{^/r{(p + s)) does not vanish for s G [0,t]. But Jacobi's 
function cn{^/r{ip + s)) vanishes at any segment of length not less than half 
of its period, thus t < T/2. By Cor. 3.3, we have T/2 < t™"^-", consequently. 
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t < t'l°'^K So the interval {0,t] does not contain conjugate points, thus the 
corresponding extremal trajectory q{s) is locally optimal (see Propos. 2.3). 

(2) Let k e (0,fco], and let the arc F contain exactly one inflection point. 
Then the function has exactly one root at the segment s G [0, t], thus t < T. 
By Cor. 3.3, we have T < tl"'^^ sot < t'^'^\ and the elastica F is locally optimal. 

(3) Let the axe F contain in its interior not less than 3 inflection points. Then 
its curvature Cg has not less than 3 roots at the interval s G (0, f). Consequently, 
the interval (0, t) contains a complete period [to, ti] of the curvature Cg such that 
Cg = at the endpoints s = to and s = ti, thus {0,t) contains a bigger segment 
with the same center: 

3[to-e,ti + s]c{0,t), s>0, 

y/F{ip + to) = K + 2Kn, y/r{ip + ti) = 5K + 2Kn, n e Z. 

So the arc F contains inside itself the elastica F = {7^ | s G [to — e,ti+ e]}. Now 
we show that the arc F is not locally optimal, this would mean that the arc F 
containing F is not locally optimal as well (indeed, if a trajectory q{s), s G [0, t], 
is locally optimal, then any its part q{s), s G [tl,t\] C [0,t] is locally optimal as 
well). 

For the arc F we have the following: 

{h + e) - {to -e)= 4K/s/r + 2e = T+2e, 

T = iiV¥{ip + to-e) + V?{ip + ti+ e))/2 = ^K + 2Kn, 

z = sn^ T = 1, Ji = ao + ai + 02, 

see (3.21). By Corollary 3.1, we have Pi""-" = min(2/r,pf ) < 2K, thus ti""-* < 

4:K/y^ = T. Consequently, (fi + e) - (to - e) = T + 2e > tl°^^ , and the interval 
(to — s,ti + e) contains a point conjugate to the instant to — e. Thus the 
arc F is not locally optimal, the more so the arc F is not locally optimal. □ 

The mathematical notion of local optimality of an extremal trajectory q[s) = 
{xs,ys,Os) w.r.t. the functional of elastic energy corresponds to stability of the 
corresponding elastica (a;,,?/,). Item (3) of Corollary 3.4 has a simple visual 
meaning: one cannot keep in hands an elastica having 3 inflection points inside 
since such an elastica is unstable. 

Remark. In the cases not considered in items (1)^(3) of Corollary 3.4, one can 
find both examples of locally optimal and non-optimal elasticae. 

Let k > ko- li z = sn^ r = 1 (i.e., the elastica is centered at its inflection 
point), then by Corollary 3.1, we have 

= min(2if,pr), G {p\,p\) C {KAK). 

For p < K we get p < the corresponding elastica contains one inflection 

point and is locally optimal, see Fig. 5. For p^^ < 2K (i.e., for k G (fc, 1), 
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Figure 5: Locally optimal elastica 
with 1 inflection point 



Figure 6; Locally non-optimal elas- 
tica with 1 inflection point 



k « 0.998) and p e {p'l^ , 2K) we get p > p™"-' = p^^ , the corresponding elastica 
contains one inflection point and is not locally optimal, see Fig. 6. 

Let k < ko and z = sn^ r = (the elastica is centered at its vertex). Then 
pConj _ pi g 3K) . If p G {K, 2K) , then p < pi, then the corresponding 
elastica is locally optimal and contains 2 inflection points, see Fig. 7. 

Let fc > fco and = sn^ r = 0, then Pi""-" = pi e {K, 2K). U p > p\), then 
p > , then the corresponding elastica is not locally optimal and contains 2 
inflection points, see Fig. 8. 




Figure 7: Locally optimal elastica Figure 8: Locally non-optimal elas- 

with 2 inflection points tica with 2 inflection points 

Corollary 3.1 provides the following description of elastica centered at inflec- 
tion points or vertices and terminating at conjugate points. 

Corollary 3.5. Let A € A^i, and let q{s) = Exp^(A), s G [0,t], be the corre- 
sponding inflectional elastica. 

(1) // the elastica q{s) is centered at its vertex (i.e., snr = 0), then the 
terminal instant t is a conjugate time iff 



(2) If the elastica q{s) is centered at its inflection point (i.e., cnr =0^, then 
the terminal instant t is a conjugate time iff 



P=^eK|neN}u{p; 



m e N}. 
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Figures 9, 10 illustrate respectively cases (1), (2) of Cor. 3.5. 



Figure 9: Conjugate point, snr = 

o,p=pr\k,T)=pi{k) 



Figure 10: Conjugate point, cnr 
0,p = pr\k,T) = 2K 



3.3 Upper bound on cut time 

On the basis of results on local optimality obtained in this section, we can 
improve the statement on upper bound on the time where elasticae lose their 
global optimahty (i.e., on the cut time tcut('^))) see Th. 12.1 [11]. The argument 
uses the obvious inequality 

icut(A) < friX), 

which follows since if a trajectory is not locally optimal, the more so it is not 
globally optimal. 

Theorem 3.3. Let A e Ni. Then tcut(A) < t(A). 

Proof. We have to prove that the extremal trajectory q{s) = Expg(A) is not 
optimal on any segment of the form s G [0,t(A) + s], e > 0. Compute the 
number r = ^{"^ip + t(A)) for the covector A = (fc, (p, r). 

Consider first the case k G (0,A;o], then t(A) = If cnr snr 7^ 0, then 
the equality icut(A) < t(A) was proved in item (1) of Th. 12.1 [11]. If cnr = 0, 
then the instant t(A) is a conjugate time by Cor. 3.5, so the trajectory q{s) is 
not locally optimal after this instant. Finally, if snr = 0, then the instant t(A) 
is a Maxwell time by item (1.1) of Th. 11.1 [11]. 
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In the case A; G (A;o, 1) wc have t(A) = and the argument is similar. If 
cnr snr 7^ 0, then the statement was proved in item (1) of Th. 12.1 [11]. If 
snr =0, then the instant t(A) is a conjugate time by Cor. 3.5. And if cnr = 0, 
then the instant t(A) is a Maxwell time by item (1.2) of Th. 11.1 [11]. □ 

4 Conjugate points on non-inflectional elasticae 

In this section we prove that inflectional elasticae (A e N2), critical elasticae 
(A € N3), and circles (A G iVe) do not contain conjugate points. 

Let A e N^- Similarly to Sec. 3, we first compute exphcitly the Jacobian of 
the exponential mapping using parametrization of extremals obtained in [11]: 

d{xt,yt,et) ^ 1 d{xt,yt,MW) ^ 32 
d{'tp,k,r) V?'cos(6it/2) d{ip,k,^) (1-/^2)^2^3/2^2 2, 



(4.1) 

J2 = C2Z^ + ciz + Co, z=snVe[0, 1], (4.2) 
p=yf?tl{2k), T = y/r{2ilj + t/k)/2, A = 1 - fc^sn^p sn^T, (4.3) 
C2 = fc^'snp cnpxi, (4.4) 



Xi = 2cnp sup E^(p) + (dnp(3 — 6sn^p) — (2 — k^)pcnp snp)E^(p) 
+ (dnp(fc^ - 2)p(l - 2sn^p) + cnp snp{k^{2p'^ - 1 + 6sn^p) 
- 2(2+p2)))E(p) + dnp{2k'^cn^psn^p + {I - k'^)p'^{2sn^ p -1)) 
+ pcnp snp (2(2 +p2) - fc^(3 + (3 - k'^)p^ + (2 - fc^) sn^p)), 

Co = -kf2{p,k) X2, (4.5) 

X2 = dnp E^(p) — fc^ cnp snp ^{p) — (1 — k'^)p'^ dnp , 

/2(p,fc) = 2(dnp(2 - fc2)p- 2E(p)) + fc^snp cnp)/k, 

C2 + ci + Co = (1 - fc2)co. (4.6) 
4.1 Preliminciry lemmas 

Lemma 4.1. For any p > and k G (0, 1), we have cq < and cq + ci +C2 < 0. 

Proof. In view of decomposition (4.5) and equality (4.6), it is enough to show 
that 

f2{p,k)>0, X2>0 Vp>0 VA;e(0,l). (4.7) 

We have 

f Mj>y\' ^ ^4 cn2p^n2p^ 
V dnp J dn^p ' 

this identity means that /2(p)/dnp increases w.r.t. variable p. But /2(0) = 0, 
thus f2{p) > for all p > and k € (0, 1). 
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Further, from the equalities 

/ X2ip) y ^ (l-fc^)(E(p)-pdn^p)2 
\dnpE{p)J dn^pF.'^ip) 

E(p) - pdn'p = ^A;V + o(/) ^ 0, 
it follows that X2{p)/{dnp E{p)) increases w.r.t. p. Then the asymptotics 

implies that X2 > for all p > and k € (0, 1). 

Inequalities (4.7) are proved, and the statement of this lemma follows. □ 

Lemma 4.2. For any n £ N, k E {0, 1), z E [0, 1] . we have J2{Kn, z, k) < 0. 

Proof. Fix any n, k, p — Kn according to the condition of this lemma. It follows 
from decomposition (4.4) that C2 = 0. Thus the function J2{z) becomes linear: 
J2{z) — ciz + Co, z E [0, 1]. By virtue of Lemma 4.1, this linear function is 
negative at the endpoints of the segment z € [0, 1]: 

J2(0) = Co < 0, J2(l) = Ci + Co = C2 + Ci + Co < 0, 

thus it is negative on the whole segment [0, 1] as well. □ 

Lemma 4.3. For any pi > there exists k = k{pi) > such that for all 

k G (0, ft), p e (0,pi), z e [0. 1] we have J2(p, z, k) < 0. 

Proof. In order to estimate the function J2 for small k, we need the correspond- 
ing asymptotics as A: — » 0: 

co = k^coo + o{k^), ci = /ci%o + o(A:'°), C2 = k^^C2o + o{k^^), (4.8) 
Coo = -Cio = -Y^(4P - sin4p) Coi(p), (4.9) 
coi = 4p^ — 1 + cos 4p + p sin 4p, 
C20 = ^^cospsinpc2i(p), 

C21 = —3 cos 2p — 48p^ cos 2p + 3 cos 6p + A2p sin 2p — 64p^ sin 2p + 2p sin 6p, 
and the asymptotics as (p, fc) (0,0): 

CO = + °(*V), C2 = -^fciy^ + o(ftiV^), (4.10) 

lib 4725 

all these asymptotic expansions are obtained via Taylor expansions of Jacobi's 
functions, see [10]. 

(1) From the equalities 

coiV 2(sm2p-2jK:os2p)^ 128 6 g 

p / p^ 45 
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it follows that coi(p) > for p > 0, whence in view of decomposition (4.9) we 
obtain that coo(p) < for all p > 0. 
Fix an arbitrary number pi > 0. 

(2) Choose any po S (0,pi). We show that there exists fcoi = A;oi(po)Pi) € 
(0,1) such that 

J2{p,z,k)<o Vpe bo,Pi] V^e [0,1] Vfce (o,fcoi). (4.ii) 

Taking into account equalities (4.8), we obtain a Taylor expansion as A; — > 0: 

J2{p,z,k) = k^Coo{p) + j^-^^{p,z,k), pe\po,Pi], ze[0,l], fce(0,fc). 
By continuity of the corresponding functions, we conclude that 
coo(p) < -m, m = m{po,Pi) > 0, 

1 J2 



10! afcio 



ip,z,k)<mi, mi = mi(po,Pi) > 0, 



whence J2 < k^{—m+ k'^nii) < for fc^ < feg^ = m/mi > 0. Inequality (4.11) 
follows. 

(3) Prom asymptotics (4.10) and equality (4.6) we conclude that 

J2 = -^kV + 0{kV), {p,k)^0. 

Thus 

3p'o>0 3k'o>0\/ p€ {0,p'o] \/k e (0, k'o) \/ z € [0, 1] J2ip, z, k) < 0. 

(4) Take j^g e (0,pi) and fcg e (0, 1) according to item (3) of this proof. Find 

^01 = fcoi(Po'Pi) according to item (2). Set fc(pi) = min(fcQ,fcoi) > 0. Then 
for any k G (0,fc(pi)) we get the following: if p € {0,Po], then J2 < by item 
(3), and if p e [po^Pi], then J2 < by item (2). So the number k{pi) satisfies 
conditions of this lemma. □ 

4.2 Absence of conjugate points 
on non-inflectional elcisticae 

Theorem 4.1. If X G N2, then the normal extremal trajectory q{t) = Exp((A) 
does not contain conjugate points for t > 0. 

Proof. In view of the symmetry i : N2 — > N2 , see [11], it is enough to consider 

the case A S -/V^. 

Denote = A. Fix any n € N and prove that the trajectory 

gi(t)=Exp,(Ai), \' = {ip,k\r)eN+, 

does not contain conjugate points t e (0,fi], tl = 2k^K{k^)nl ^/r. 
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Consider the family of extremal trajectories 

Y = {q\t)=E^p,iX^)\te[0,t{]}, 

X' = {ip,k',r) G N+, tl = 2k'K{k')n/Vr, se[0,l], 

where the covector = {(f,k^,r) is equal to A given in formulation of this 
theorem, and covector A° = {(p. fc". r) will be chosen below in such a way that 
the parameter k^ is sufficiently small. 

According to the Lemma 4.3, choose a number k{p^) G (0, 1) correspond- 
ing to the number = K{k^)n. Choose any e {0,k{p^)) and set A° = 
{ip,k°,r)GN+. 

By Lemma 4.3, for any p G {0,p^] and any z G [0, 1] we have J2{p, z, /c") < 0. 
By Lemma 4.2, for any z G [0, 1] and any k G [ko, ki] we have J2{K{k)n, z, k) < 
0. 

Taking into account equahty (4.1) and relations (4.2), (4.3), we conclude 
that the trajectory 7" docs not have conjugate points at the segment t € {0,ti], 
and for any trajectory 7", s G [0, 1], the endpoint t = tl is not conjugate. Now 
the statement of this theorem follows from Corollary 2.2. □ 

4.3 Absence of conjugate points for special cases 

The absence of conjugate points on for extremals Xt G N2 implies a similar fact 

for XtGNsU Nq. 

Theorem 4.2. If X G N^U Nq, then the extremal trajectory q{t) = Exp((A) 
does not contain conjugate points for t > 0. 

Proof. Let A e A^3 U A^e • Since the set A^3 U Nq belongs to the boundary of the 
domain N2, one can construct a continuous curve A* : [0, 1] — > AT such that 
A* G N2 for s G [0, 1) and A^ = A. 

Consider the family of extremal trajectories q^{t) — Expj(A*), t > 0, s G 
[0,1]. It follows from Th. 4.1 that for s G [0,1) the trajectory (/"(t) docs not 
contain conjugate points t > 0. Then we conclude from Corollary 2.3 that the 
trajectory q^{t) = Expj(A) does not contain conjugate points for t > 0. □ 

5 Final remarks 

We sum up the study of conjugate points in Euler's elastic problem. 

Any inflectional elastica contains an infinite sequence of isolated conjugate 
points. Visually, the first conjugate point occurs between the first and third 
inflection points. More intrinsically, the flrst conjugate point belongs to the 
interval (T/2,3T/2), where T is the period of oscillation of the pendulum that 
parametrizes the vertical subsystem of the normal Hamiltonian system. The 
first conjugate point is contained in the arc bounded by the first Maxwell points 
corresponding to discrete symmetries [11]. 
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Non-inflcctional and critical elasticae, circles and straight lines do not con- 
tain conjugate points. 

On the basis of this information about conjugate points and the description 

of Maxwell points obtained in [11]. one can study the global structure of the 
exponential mapping in Euler's elastic problem; describe the domains where the 
exponential mapping is diffeomorphic, and find a precise characterization of cut 
points. Another interesting question for further study is the structure of the 
caustic in Euler's problem. This will be the subject of our forthcoming work. 
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